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LS 1+61 303: the hadronic contribution revisited 
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Abstract LS 1+61 303 has been detected by the Che- 
renkov telescope MAGIC at very high energies, present- 
ing a variable flux along the orbital motion with a max- 
imum clearly separated from the periastron passage. In 
the light of the new observational constraints, we revisit 
the discussion of the production of high-energy gamma 
rays from particle interactions in the inner jet of this sys- 
tem. The hadronic contribution could represent a major 
fraction of the TeV emission detected from this source. 
The spectral energy distribution resulting from pp inter- 
actions is recalculated. Opacity effects introduced by the 
photon fields of the primary star and the stellar decretion 
disk are shown to be essential in shaping the high-energy 
gamma-ray light curve at energies close to 200 GeV. We 
also present results of Monte Carlo simulations of the 
electromagnetic cascades developed very close to the pe- 
riastron passage. We conclude that a hadronic micro- 
quasar model for the gamma-ray emission in LS I +61 
303 can reproduce the main features of its observed high- 
energy 7-ray flux. 

Keywords X-ray binaries • mass loss and stellar winds • 
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1 Introduction 

LS I +61 303 is a Be X-ray binary with strong and 
variable radio emission. It has been suggested as the 
potential counterpart of the gamma-ray sources 2CG 
135+01 and 3EG J0241+6103 (Gregory & Talor 1978, 
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Kniffen et al. 1997). Massi et al. (2001) detected rela- 
tivistic jets in the source, which was then included in the 
microquasar class. The existence of jets was confirmed 
by the observations presented by Massi et al. (2004). 
Very recently, LS I +61 303 was detected by MAGIC, a 
very large atmospheric imaging Cherenkov telescope lo- 
cated at La Palma (Albert et al. 2006). This makes LS I 
+61 303 the second high-energy microquasar detected by 
ground-based telescopes. Moreover, MAGIC found clear 
evidence for variability in the detected signal: the source 
was stronger at orbital phases 0.5-0.6, whereas the peri- 
astron passage occurs at phase 0.23. 

At X-rays the source has been detected by XMM- 
Newton, Beppo-SaX and INTEGRAL (Sidoli et al. 2006, 
Chernyakova et al. 2006) showing variability on short 
timescales and a harder spectrum when the source was 
brighter (Sidoli et al. 2006). 

The orbital parameters of LS I +61 303 have been 
determined by Casares et al. (2005). The eccentricity of 
the system is 0.72+0.15 and the orbital inclination is 
i ~ 30 ± 20 deg. The nature of the compact object is not 
well-established. A low-mass black hole cannot be ruled 
out (Massi 2004). If a mass of ~ 2.5 M Q is assumed for 
the compact object, then the mass of the Be star should 
be ~ 12 M Q , for i = 30 deg. 

From a theoretical point of view, both accreting and 
non-accreting models have been proposed to explain the 
high-energy emission of LS 1+61 303. In the accreting 
scenario (e.g. Bosch- Ramon et al. 2006a) a relativistic jet 
is launched from the surroundings of the compact object. 
Strong magnetic fields are supported by an underlumi- 
nous accretion disk, which is advection-dominated. High- 
energy emission can be produced in the jet by leptonic 
(e.g. Bosch-Ramon et al. 2006b, Gupta & Bottcher 2006, 
Bednarek 2006a) or hadronic (Romero et al. 2003, 2005) 
interactions. In non-accreting models (Dubus 2006b, Cher- 
nyakova et al. 2006) particles are accelerated up to rel- 
ativistic energies in the interacting region where a pul- 
sar wind collides with the stellar wind. Inverse Compton 
(IC) cooling of electrons (Maraschi & Treves 1981, Dubus 
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2006b) or pp interactions (Chernyakova et al. 2006) could 
result in the production of high-energy photons. 

A colliding wind model has been proposed as the 
mechanism producing the high-energy emission in the 7- 
ray binary PSR B1259-63, a system that contains a radio 
pulsar in an eccentric orbit around a Be star (Aharonian 
et al. 2005). In the case of LS 1+61 303, the existence 
of jets extending up to ~ 400 AU from the core seems 
to support an accretion/jet model (sec the discussion in 
Mirabel 2006). Albert et al. (2006) favor a microquasar 
model and suggest that the detection of the high-energy 
signal after the periastron passage would support a lep- 
tonic origin for the emission. However, 7-ray propagation 
effects alone could be responsible for such a behavior 
(Bednarek 2006a). 

In this paper we revisit the hadronic model for high- 
energy emission in LS 1+61 303 proposed by Romero et 
al. (2005) before the MAGIC detection. Using an im- 
proved accretion model and more sophisticated calcula- 
tions for the opacity effects inside the binary system we 
show that a hadronic origin for the emission detected by 
MAGIC can not be ruled out. 



2 The model 

We consider gamma-ray production in hadronic interac- 
tions of relativistic protons from a jet launched close to 
the compact object and cold protons from the equato- 
rial wind of the Be star (Romero et al. 2003, 2005). The 
star has a slow equatorial wind that fits a density profile 
p w (r) = p (r/R*)- n , with n = 3.2 and p = 5 x 10" 11 
g cm -3 (Marti and Paredes 1995). The wind remains 
mainly near to the equatorial plane, confined in a disk 
with half-opening angle ip = 15°, and an initial outflow- 
ing radial velocity v r o ~ 5 km s~ x (Marti and Paredes 
1995). The disk effective temperature close to the star is 
2disk = 17000 K. The values of the adopted parameters 
in our model are listed in Tabled 

The pp interactions can occur either because there is 
some mixing of the stellar wind with the jet matter or 
because some protons escape from the jet into the wind. 
The level of interaction is quantified through a "mixing 
factor", here assumed as / m oc v~j ~ 0.1. This phe- 
nomenological prescription accounts for a more efficient 
rejection of particles when a larger relative velocity is 
thought to enhance the boundary effects between the jet 
and the wind 1 . At phase 4> = 0.5 the mixing factor is 
maximum, reaching f m ~ 0.5. 

The jet is assumed to be in a steady state. We are not 
dealing here with the details of the launching mechanism, 
which is supposed to be related to magneto-centrifugal 
effects (Blandford & Payne 1982). The hadronic jet power 
is a fraction of the accretion power: Lj e t = 9jet-MaccrC 2 . 

The problem of matter exchange through the boundary 
layers of a relativistic jet is a difficult one. Its proper treat- 
ment is far beyond the scope of this work. 



The value of q- 3Ct is taken as a free parameter. In order to 
reproduce the MAGIC observations we have that qj C t ~ 
0.1. The mass accretion rate M acci is strongly dependent 
on the relative velocity of the compact object with re- 
spect to the wind. To estimate the evolution of the accre- 
tion rate onto the black hole we have taken into account 
the azimuthal wind velocity ~ 1.1 x 10 7 (i?*/r)/ sini 
cm s _1 (Gregory & Neish 2002, Casares et al. 2005) in 
addition to the radial velocity, v T = 3 x 10 5 (r/i?*)"~ 2 
cm s _1 . The Bondi-Hoyle accretion regime 2 was consid- 
ered to obtain the values of M accr shown in Figure 1. 
This is just a rough approximation used by several au- 
thors (Marti and Paredes 1995, Gregory & Neish 2002). 
Close to the compact object there should be a transi- 
tion to disk accretion, but the absence of disk signatures 
in the X-ray spectrum suggests that the accretion disk 
in LS 1+61 303 is small. We shall assume that changes 
in M accr are propagated to the jet on timescales much 
shorter than the orbital period. 
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Fig. 1 Normalized mass accretion rate. At the periastron 
passage (phase 0.23) it is M acc = 1.7 x 10 19 g s _1 . 



3 The very high-energy spectrum and light 
curve 

In Figure[21we show the gamma-ray luminosity expected 
from pp — > pp + 7T° and the subsequent neutral pion de- 
cays. The light curve is plotted for photon energies ~ 200 
GeV, which is the lower energy reported by MAGIC for 
LS I +61 303. The differential 7-ray emissivity is cal- 
culated applying the <5-function approximation (Aharo- 
nian & Atoyan 2000). The jet is assumed to expand in 
a conical way (see parameters in Table The mag- 
netic field (B < 10 s G at the base of the jet) results 
from the equipartition condition (see Bosch Ramon et al. 
2006a) and then decreases with the distance as oc z~ 2 . 
Notice that the field will change with the accretion rate 
and hence, with the orbital phase. Since hadronic energy 
losses are negligible, the maximum proton energy has 
been obtained at the jet formation point equating the 

2 We notice that the close approach of the stars during the 
periastron might induce a transient Roche-lobe overflow. 
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Table 1 Model parameters 



Parameter: description [units] 


values 


Mi,: stellar mass [Mq] 


12 


7?*: stellar radius [Rq] 


10 


-Mbh: compact object mass [Mq] 


2.5 


e: eccentricity 


0.72 


i: orbital plane inclination [°] 


30 


ui: angle of periastron [°] 


20 


Rq: initial jet radius [i?sch] 


5 


zq: jet initial point in the compact object RF [-Rs c h] 


50 


X: jet semi-opening angle tangent 


0.1 


_B eq : equipartition magnetic field at the base of the jet [G] 


10 s 


a: relativistic proton power-law index 


2.5 


7™ ln : lowest Lorentz factor of relativistic protons 


2 


7jf ax : highest Lorentz factor of relativistic protons 


~ 10 7 


r: macroscopic jet Lorentz factor 


1.25 


T*: stellar surface temperature [K] 


26000 


Tdi s k: circumstellar disk temperature [K] 


17000 


^?disk: circumstellar disk outer radius [R*] 


12 


Ri n : radius of the circumstellar disk brightest region [R*] 


3 


Ldisk: circumstellar disk luminosity [erg s" 1 ] 
(p: disk half-opening angle [°] 


2 x 10 37 


15 


« r o: radial velocity of the wind at the base [cm s _1 ] 


3 x 10 5 


v^y. azimuthal velocity of the wind at the base [cm s _1 ] 


1.13 x 10 7 /sini 


Po : density of wind at the base [g cm -3 ] 


5 x 10" 11 



shock acceleration rate to the proton synchrotron rate. 
This leads to j p = 7 x 10 6 at the base of the jet during 
the periastron passage. Photohadron interactions in the 
stellar field (the typical energy of the stellar photons is 
~ 1 eV) are well below the threshold even for the most 
energetic protons. These type of losses are possible with 
X-ray photons form the accretion disk, but the X-ray 
luminosity is low enough as to neglect them in a first 
approximation. The differential distribution of relativis- 
tic protons has an index a — 2.5 in order to match the 
slope of the MAGIC spectrum. We note that the usual 
particle acceleration through the Fermi mechanism can 
lead to such a soft distribution for shocks with low Mach 
number. Actually, a could be time dependent, affected 
at some extent by the strong variations in the accretion 
rate. 

There are two maxima in the light curve shown in 
Figure [3 one at the periastron and the other around 
phase 0.5, when the relative velocity is minimum. The 
solid line shows the luminosity corrected by 77 absorp- 
tion in the stellar and disk radiation fields. We have cal- 
culated the opacity of the anisotropic stellar field as in 
Dubus (2006a), taking into account angular effects and 
the finite size of the star. The orientation of the orbit 
is given by Casares et al. (2005). The circumstellar disk 
is modeled as a blackbody of Tdisk = 17000 K for the 
inner region (up to an inner radius of 3 stellar radii) and 
then with an emissivity that goes as oc (Waters 1986, 
Bosch- Ramon et al. 2006b). The disk farther truncates 
at 12 i?*. The inner region of the disk is normalized to 
emit Ldisk ~ 2 x 10 37 erg s -1 , which is a non negligi- 
ble fraction of the total thermal emission of the system 



(Casares et al. 2005). The total 77 cross-section has been 
considered for interactions with photons from the disk. 
The energy dependence of the optical depth from both 
disk and stellar contribution is shown in Figure|21 Figure 
01 presents the variation of the total optical depth with 
the orbital phase for a photon of energy of 200 GeV. 

The spectral energy distribution (SED) frompp inter- 
actions, estimated for different orbital phases, is shown 
in Figure Close to the periastron the absorption re- 
sults to be significant, and mainly dominated by the cir- 
cumstellar disk emission. On the other hand, at phase 
(f> = 0.53 most of the radiation escapes from the source, 
making it detectable at this phase and not during the 
first accretion peak, when the ambient wind density is 
maximum. Notice that the more realistic modeling of 
the wind and the absorption changes these results from 
those previously presented by Romero et al. (2005) 3 . 

We have computed the high-energy 7-ray spectra re- 
sulting from cascades traversing the anisotropic stellar 
radiation field during the periastron passage. The star is 
characterized by a surface temperature T* = 26000 K, 
and a stellar radius i?* = 10 Rq. Monte Carlo simula- 
tions were performed after developing a computational 
code based on the scheme outlined by Protheroe (1986) 
and Protheroe et al. (1992). The 7-ray spectra produced 
through IC interactions were approximated as in Jones 
(1968), modified in similar way to Bednarek (1997). We 
introduced the effects of the finite size of the star and 



3 Also, we considered in the present case the usual conven- 
tion for the orbital phase as time proportional. In Romero 
et al. (2005) we have used the orbital parameter as linearly 
related to the so-called true anomaly. 
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Fig. 2 Calculated 7-ray luminosity from hadronic interac- 
tions. The dashed line corresponds to the generated luminos- 
ity. In solid line, we show the luminosity corrected by ab- 
sorption in the stellar and disk photon fields, according to 
the values of r shown in Figure 0] 



the spatial variation of the ambient photon field density 
by considering the geometric configuration as in Dubus 
(2006a). The magnetic field in the cascade region, orig- 
inated in the star, was assumed to be lower than 0.1 
G, assuring in this way a strongly Compton dominated 
regime. We followed the cascades induced by the injected 
gamma-ray flux with mentioned photon index a = 2.5, 
during the periastron passage (i.e. orbital separation of 
- 2.56 i?*). 

In Figure|H|we present the results of the simulations. 
The one-dimensional cascade development is calculated 
in the anisotropic stellar field along the line of sight, 
forming a 30 deg angle with the jet axis. At phase 0.5 the 
opacity is not enough to sustain a local cascade (gamma- 
rays that escape in the direction of the star can produce 
cascades in the stellar photosphere, and some gamma- 
rays could, in principle, be redirected toward the ob- 
server, see Bednarek 2006a and 2006b). We notice that 
at the periastron passage the effect of the cascades is to 
produce a softer 7-ray spectrum with index a ~ 2.8 in 
the MAGIC energy range. Such a feature should be de- 
tectable through larger exposures than those reported in 
Albert et al. (2006). 



4 Discussion 

The luminosity obtained at phase 0.5 implies that the 
constant coupling the jet kinetic power and the accretion 
luminosity should be <7j c t ~ 0.1 in order to explain the 7x 
10 33 erg s _1 inferred from MAGIC detections. Powerful 
jets can be present in microquasars, as showed by Gallo 
et al. (2005) for the case of Cygnus X-l. 

At energies below the break of the hadronic spectrum 
(~ GeV) the emission of leptonic origin should change 
the shape of the SED yielding the EGRET spectrum. 
The hadronic interactions that produce neutral pions 
also generate relativistic electrons and positrons by the 
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Fig. 3 Optical depth as function of the 7-ray photon energy. 
It was calculated for fixed values of the orbital phase (f> indi- 
cated in the figure. The lower panel shows absorption by the 
stellar photon field, and the upper one, by the photon field 
of the circumstellar disk. 



decay of charged it- mesons. We have not computed here 
the emission of such leptons, for we are concerned only 
with the very high energy range. Some estimations of the 
synchrotron luminosity of secondaries were presented in 
Romero et al. (2005). 

Both leptonic and hadronic microquasar models for 
the gamma-ray emission in LS I +61 303 can reproduce 
the main features of the observed gamma-ray emission. 
The detection of the source at phase 0.5 is a consequence 
of the opacity effects and is independent of the dominant 
radiative process in the jet. Even in hadronic models 
like the one presented here, near the periastron passage 
gamma-rays will initiate electromagnetic cascades with 
IC emission, hence the actual situation could be a mix- 
ture of radiative processes. More realistic simulations of 
the photospheric cascades should include the effects of 
circumstellar disk absorption. 

In a pure leptonic model the radiative losses equating 
the acceleration rate fix the maximum energies of the 7- 
ray photons at ~ 1 TeV (Bosch- Ramon et al. 2006b). 
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Fig. 4 Optical depth for the propagation of 7-ray photons 
with energy E = 200 GeV. 
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Fig. 5 Estimated spectral energy distribution, affected by 
absorption, for different values of the orbital phase. In this 
plot qjet = 1. Note that at tj> = 0.53 (solid line) the expected 
signal should be stronger than at the periastron passage (</> = 
0.23) for the range of energies covered by MAGIC. The actual 
cutoff for high-energy photons is around 1 PeV. 
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Fig. 6 Emerging spectral energy distribution after the IC 
dominated cascade developed in the stellar photon field (the 
effect of the decretion disk of the Be star is not included). The 
straight line only indicates the slope of the injected spectrum, 
with photon index a = 2.5 (which mimics the original proton 
spectral index). These results were obtained through Monte 
Carlo simulations, properly taking into account the geometric 
configuration. The curves were normalized with gj e t = 1 and 
the actual cutoff is around 1 PeV. 

is that the high accretion rates and the observed low X- 
ray luminosity seem to exclude a scenario based upon 
an accreting pulsar, since there is no trace of the emis- 
sion from the heated surface. Future observations with 
MAGIC will help to detect the source close to the perias- 
tron and the spectral evolution along the orbit, providing 
in this way more constraints to the models. 
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Detection of 7-ray emission of LS I +61 303 at higher 
energies than those observed by MAGIC (e.g. at E > 10 
TeV) could give support to a hadronic model, and to the 
presence of intense magnetic fields at the base of the jet. 

A clear distinction about the nature of most of the 
primary gamma-rays can be established through the de- 
tection of neutrinos, which are produced only in hadronic 
models. In the context of our model, we expect around 
3-4 muon neutrinos per squared km per year at Earth, 
so the source might be detectable by IceCube (Chris- 
tiansen et al. 2006, see also the discussion in Torres & 
Halzen 2006). 

As we mentioned in the Introduction, an alterna- 
tive model for the gamma-ray production in LS I +61 
303 involves an energetic pulsar which could generate 
a strong wind that would stop the accretion. Particles 
might be accelerated at the colliding wind region, pro- 
ducing gamma-rays through IC and pp interactions (Dubus 
2006b, Chernyakova et al. 2006). However, it is not clear 
how the morphological and variability properties of the 
system might be explained in this scenario. What is clear 
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